Precise determination of the /o(600) and /o(980) pole parameters from a dispersive 

data analysis 
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We use our latest dispersive analysis of tttt scattering data and the very recent Ke4, experimental 
results to obtain the mass, width and couplings of the two lightest scalar- isoscalar resonances. These 
parameters are defined from their associated poles in the complex plane. The analytic continuation 
to the complex plane is made in a model independent way by means of once and twice subtracted 
dispersion relations for the partial waves, without any other theoretical assumption. We find the 
/o(600) pole at (457li3) " ^(2791:^) MeV and that of the /o(980) at (996 ± 7) - 1(2516") MeV, 
whereas their respective couplings to two pions are 3.59lo'i3 GeV and 2.3 ± 0.2 GeV. 

PACS numbers: 14.40.Be,11.55.Fv,13.75.Lb 



The /o(600) or sigma and /o(980) resonances are of 
great interest in several fields of Physics. First, the two 
pion exchange in the scalar-isoscalar channel, 1=0, J=0, 
where these resonances appear, plays a key role in Nu- 
clear Physics, where the nucleon-nucleon attractive in- 
teraction has been for long [1] modeled by the exchange 
of a "sigma" resonance. Second, this channel is also rel- 
evant for the QCD non-Abelian nature, since it is where 
the lightest glueball is expected to appear. However, 
the glueball identification is complicated by its possi- 
ble mixing into different states, like the /o(600), /o(980) 
and heavier /o resonances, which may be qq mesons, 
tetraquarks, molecules, or most likely a mixture of them 
all. Actually, most of the controversy around these reso- 
nances comes from the identification of scalar multiplets 
— see the Review of Particle Physics (PDG) "Note on 
Scalar Mesons" H. Third, the /o(600), being the hght- 
est hadronic resonance with vacuum quantum numbers, 
plays a relevant role in many models of QCD sponta- 
neous chiral symmetry breaking. Furthermore, this state 
is of interest in order to understand why, despite being 
so light and strongly coupled to pions, it plays such a 
small role, if any, in the saturation [3j of the low energy 
constants of Chiral Perturbation Theory (ChPT). More- 
over, the position of this pole could be setting the limit of 
applicability of the chiral expansion. Finally, this state is 
of interest for Electroweak Physics due to its many sim- 
ilarities — but even more by its many differences — with 
the Higgs mechanism now under scrutiny at LHC. 

Still, the properties of these resonances are the sub- 
ject of an intense debate. Let us recall that the a was 
listed in the PDG as "not well established" until 1974, 
removed in 1976, and listed back in 1996. This was due 
to its width being comparable to its mass, so that it 
barely propagates and becomes a broad enhancement in 
the traditional, and often contradictory, tttt scattering 
analyses, extracted from TriV — > ttttA'^ experiments, using 
different models affected by large systematic uncertain- 



ties. After 2000 these resonances have been observed in 
decays of heavier mesons, with well defined initial states 
and very different systematics from tttt scattering, which 
led the PDG to consider, in 2002, the /o(600) as "weU 
established" , but keeping until today a too conservative 
estimate of: "Mass: 400 to 1200 MeV" and "Width: 600 
to 1000 MeV". For the /o(980) the situation is not much 
better, with an estimated width "from 40 to 100 MeV" . 
However, not all the uncertainty comes from experiment. 
The shape of these resonances varies from process to pro- 
cess and that is why their masses and widths are quoted 
from their process independent pole positions, defined as 
^Spoic ^ M — iT/2. But many models do not implement 
rigorous analytic continuations and lead to incorrect de- 
terminations when poles are deep in the complex plane 
or close to threshold cuts, as it happens with the /o(600) 
and the /o(980), respectively. Actually, this is one of the 
main causes of the huge PDG uncertainties [4I. 

This model dependence can be avoided by using dis- 
persive techniques, which follow from causality and cross- 
ing, and provide integral relations and a rigorous analytic 
continuation of the amplitude in terms of its imaginary 
part in the physical region, which can be obtained from 
data. For example, dispersion relations combined with 
ChPT determine the a polecat 440 - i 245 MeV Q or 
(470±50)-i (260±25) MeV iSij. We focus here on disper- 
sive analyses, but othe r ap proaches yield similar values 
[SIlZI — see Table H] and [14| for a review and references. 

Generically, the main difficulty lies in the calculation 
of the left cut integral, which in [J, |5| was just approx- 
imated. This left cut is due to crossing symmetry and 
can be incorporated rigorously in a set of infinite cou- 
pled equations written long ago by Roy |15j ( see also 
[16| for applications and references). Recently, Roy equa- 
tions have been used to study low energy tttt scattering 
jlTj . sometimes combined with ChPT [18], or also to test 
ChPT [19| , as well as to solve old data ambiguities [20J . 
Most recently [§], the /o(600) and /o(980) poles were 
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TABLE I: Other recent determinations of tlie a pole and cou- 
pling, using analyticity. Results come from Roy eqs. and 
ChPT 8], conformal fits to K14 decays and averaged tttt data 
around 800-900 MeV with only statistical [3] or also system- 
atic [101 uncertainties, the chiral unitary approach [ij] ( only 
statistical error), a iC-matrix with a form factor shap e [l^ . 
and ChPT±elastic dispersion relations (two- loops [13l]). 
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TABLE IL Recent determinations of /o(980) parameters. For 
[22 ] our estimate covers the six models considered there. The 
last three poles come from scattering matrices and the rest 
from production experiments. 



shown to lie within the apphcability region of Roy eqs. 
Since data were not rehable and to improve accuracy, 
Roy eqs. were supplemented by ChPT predictions in [Sj, 



to yield : ^ = 441Ig - i 272119 5 MeV, without using 



data below 800 MeV on S and P waves. In that work 
an /o(980) pole is also found at ^i = 1001 - i 14 MeV. 
Note that, generically, tttt scattering data around 900 
MeV tend to produce a narrower /o(980) 0,11, EH than 
that seen in production processes or the PDG estimate. 
In Table |TT] we list some other recent determinations of 
the /o(980) parameters. 

Our aim in this work is to provide a precise and model 
independent simultaneous determination of the /o(600) 
and /o(980) parameters from data alone, profiting from 
two relevant results developed over the last half year: 
On the one hand, the final analysis of Kgi decays by the 
NA48/2 Collaboration [26|, which provides reliable and 
precise tttt scattering phases below the mass of the kaon. 
On the other hand, a set of Roy-like eqs. — called GKPY 
eqs. and developed by our group [27[ — which is much 
more stringent in the resonant region than standard Roy 
eqs. The reason is that, in order to avoid divergences, dis- 
persion relations are weighted at low energy with "sub- 
tractions" , but then amplitudes are only determined up 
to a polynomial, whose coefficients depend on thresh- 
old parameters. Since Roy eqs. have two subtractions 
they have an s polynomial term multiplied by the isospin 



2 scalar scattering length, whose large uncertainty thus 
grows markedly in the /o(600) and /o(980) region. In 
contrast, the GKPY eqs. have just one subtraction and 
their output, even without using ChPT predictions at all, 
provides [27[ a very precise description of tttt scattering 
data, discarding a long-standing conflict concerning the 
inelasticity — and to a lesser extent the phase shift — right 
above the /o(980) region. 

If we now use these GKPY dispersion relations to con- 
tinue analytically that amplitude, we find: 
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(457^13) - «(279i") MeV 
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Let us describe next the whole approach in detail and 
provide determinations for other quantities of interest, 
like their couplings and the p{770) parameters, as well as 
other checks of our calculations from Roy eqs. 

Ours is what is traditionally called an "energy- 
dependent" analysis of tttt scattering and Kg4 decay 
data |28l. | 2 9l| — in particular the latest results from 
NA48/2 2611 ■ Our procedure, described in a series of 
works [271 130| was first to obtain a simple set of uncon- 
strained fits to these data (UFD) for each partial wave 
separately up to 1420 MeV, and Regge fits above that 
energy. Next we obtained constrained fits to data (CFD) 
by varying the UFD parameters in order to satisfy within 
uncertainties two crossing sum rules, a complete set of 
Forward Dispersion Relations as well as Roy and GKPY 
eqs., while simultaneously describing data. The details 
for all CFD waves can be found in [27| , but since we are 
now interested in the scalar isoscalar partial wave ig , we 
show in Fig.[T]the resulting 6q phase shift. It should be 
noticed that the CFD result is indistinguishable to the 
eye from the UFD, except in the 900 to 1000 MeV region, 
which we also show in detail and is essential for the de- 
termination of the /o(980) parameters. Note that both 
the UFD and CFD describe the data in that region, but 
the GKPY dispersion relations require the CFD phase 
to lie somewhat higher than the UFD one. This is rele- 
vant since it yields a wider /o(980), correcting the above 
mentioned tendency to obtain a too narrow /o(980) from 
unconstrained fits to tttt scattering data alone. In the in- 
ner top panel, we show the good description of the latest 
NA48/2 data on Ke4 decays, which are responsible for 
the small uncertainties in our input parametrization and 
constrain our subtraction constants. As seen in Fig. [l] 
the inelasticity tjq ' shows a "dip" structure above 1 GeV 
required by the GKPY eqs. [27|, which disfavors the al- 
ternative "non-dip" solution. Having this long-standing 
"dip" versus "no-dip" controversy [3l| settled [27| is very 
relevant for a precise /o(980) determination. 

The interest of this CFD parametrization is that, while 
describing the data, it satisfies within uncertainties Roy 
and GKPY relations up to their applicability range, 
namely 1100 MeV, which includes the /o(980) region. 
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FIG. 1: SO wave phase and inelasticity from UFD and CFD. 
Dark bands cover the uncertainties. Data come from [26i. i28||. 



In addition, the three Forward Dispersion Relations are 
satisfied up to 1420 MeV. In Fig. [5] we show the fulfill- 
ment of the SO wave Roy and GKPY eqs. and how, as 
explained above, the uncertainty in the Roy eq. is much 
larger than for the GKPY eq. in the resonance region. 
The latter will allow us now to obtain a precise determi- 
nation of the /o(600) and /o(980) poles from data alone, 
i.e. without using ChPT predictions. 

Hence, we now feed our CFD parameterizations as in- 
put for the GKPY and Roy eqs., which provide a model 
independent analytic continuation to the complex plane, 
and determine the position and residues of the second 
Riemann sheet poles. It has been shown [8| that the 
/o(600) and /o(980) poles lie well within the domain of 
validity of Roy equations, given by the constraint that the 
t values which are integrated to obtain the partial wave 
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FIG. 2: Fulfillment of SO wave Roy and GKPY equations. 
The CFD parametrization is the input to both the Roy and 
GKPY eqs., and is in remarkable agreement with their output. 
Note how the uncertainty in the Roy eq. is much larger than 
that of the GKPY eq. above roughly 500 MeV. 



representation at a given s should be contained within a 
Lehmann-Martin ellipse. These are conditions on the an- 
alytic extension of the partial wave expansion, unrelated 
to the number of subtractions in the dispersion relation, 
and equally apply to GKPY eqs. 

Thus, in Table Hn] we show the /o(600), /o(980) 
and p(770) poles resulting from the use of the CFD 
parametrization inside Roy or GKPY eqs. We consider 
that our best results are those coming from GKPY eqs. 
since their uncertainties are smaller, although, of course, 
both results are compatible. 

Several remarks are in order. First, statistical un- 
certainties are calculated using a MonteCarlo Gaussian 
sampling of the CFD parameters with 7000 samples dis- 
tributed within 3 standard deviations. A systematic 
uncertainty due to the different charged and neutral 
kaon masses is relevant for the /o(980) due to the ex- 
istence of two KK thresholds separated by roughly 8 
MeV, which we have treated as a single KK threshold 
at iriK = (w^o — mx+)/'^ — 992 MeV. In order to es- 
timate this systematic uncertainty, we have refitted the 
UFD and CFD sets to the extreme cases of using ttiko or 
TOj^+ instead of m^ . As it could be expected, the only 
significant variation is for the /o(980) — actually, only for 
its half width, which changes by ±4.4 MeV for GKPY 
eqs, and ±5.6 MeV for Roy eqs. The /o(600) changes 
by roughly 1 MeV and the /9(770) barely notices the 
change — less than 0.1 MeV. The effect on residues is 
smaller than that of rounding the numbers. We have 
added all these uncertainties in quadrature to the statisti- 
cal ones. Second, both the mass and width of the /o(600) 
are compatible with those in ref. '8* within one standard 
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TABLE III: Poles and residues from Roy and GKPY eqs. 



deviation. Since we are not using ChPT and ref. 8] did 
not use data below 800 MeV, this is a remarkable check 
of the agreement between ChPT and low energy data. 
Third, the /o(980) width is no longer so narrow — as it 
happens in typical tttt scattering analyses — and we find 
r = 50+22 MeV, very compatible with results from pro- 
duction processes. The mass overlaps within one stan- 
dard deviation with the PDG estimate. These results 
show that the effect of the too narrow /o (980) pole and 
the use of further theoretical input like ChPT do not 
affect significantly the resulting /o(600) parameters. 

In Table Hill we also provide for each resonance its cou- 
pling to two pions, defined from its pole residue as: 



g = — IGtt lini (s 



Spoie)U{s)i2£+l)/i2py' (3) 



where p^ = s/4 — mj. This residue is relevant for mod- 
els of the spectroscopic nature of these particles, particu- 
larly for the /o (600) [32| , which are beyond the pure data 
analysis scope of this work. Differences between previous 
values of these couplings can be seen in Tables H] and HT] 
In summary, using a recently developed dispersive for- 
malism, which is especially accurate in the resonance re- 
gion, we have been able to determine, in a model inde- 
pendent way, the /o(600), /o(980) poles and couplings 
from data with no further theoretical input. We hope 
this works helps clarifying the somewhat controversial 
situation regarding the parameters of these resonances. 
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